Study of light-assisted collisions between a few cold atoms in a microscopic dipole trap 
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We study light-assisted collisions in an ensemble containing a small number (~ 3) of cold ^'^Rb 
atoms trapped in a microscopic dipole trap. Using our ability to operate with one atom exactly 
in the trap, we measure the one-body heating rate associated to a near-resonant laser excitation, 
and we use this measurement to extract the two-body loss rate associated to light-assisted collisions 
when a few atoms are present in the trap. Our measurements indicate that the two-body loss rate 
can reach surprisingly large values /? > 10~*cm"^s~^ and varies rapidly with the trap depth and the 
parameters of the excitation light. 

PACS numbers: 34.50.Cx,34.50.Rk,37.10.Gh 
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I. INTRODUCTION 

Extensive experimental and theoretical studies have 
been devoted in the last decades to light-assisted colli- 
sions, using cold atoms held in a magneto-optical trap [l|, 
including at the few atom level [2|, or in large optical 
dipole traps [3|, |J| . In small dipole traps with size com- 
parable to the wavelength of the light, light-assisted col- 
lisions are used to prepare or probe mesoscopic atomic 
ensembles, opening new avenues in condensed matter 
physics and quantum information processing. For ins- 
tance, they are at the heart of the preparation of in- 
dividual atoms in microscopic optical dipole traps \5r 
0], standing waves Q or three dimensional optical lat- 
tices [9|. They are also at the origin of sub-poissonian 
atom number distributions in a mesoscopic atomic en- 
semble [10| . Finally, they constitute an important tool to 
understand quantum phases, as demonstrated recently 
with atoms in optical lattices UJ, |l2|. While concep- 



tually simple, the theoretical description of light-induced 
collisions is known to be cumbersome due to the com- 
plex interplay between atomic multi-level structure and 
atom-light coupling. The situation is even worse when 
considering tightly-confined atomic ensembles where the 
trapping potential acts on the same length scale as the 
interaction between the atoms. As a consequence no theo- 
retical prediction for the loss rates and their dependency 
on parameters such as the atomic density or the light pa- 
rameters is available to date for this system. The absence 
of reported measurements makes the situation even more 
dramatic. 

In this paper, we report on an experimental study 
of light-assisted collisions between cold atoms that are 
tightly confined in a microscopic dipole trap. To allow 
for future theoretical modeling of our data, we imple- 
mented as closely as possible the gedanken experiment 
where merely two atoms in the ground state (here, *^Rb 
in state 551/2 F = 2) collide in the presence of a nearly- 
resonant laser field. In our case, the loading of the dipole 
trap is non deterministic [10| and we operate with a ty- 



pical average atom number of ^ 3. We then illuminate 
the trapped atoms with a pulse of near-resonant light 
with known frequency and intensity, in order to trigger 
losses. The near-resonant light has two effects : it heats 
the atoms individually out of the trap and it induces two- 
body losses, which we wish to study. To separate the two 
contributions, we proceed in two steps. First, we use our 
ability to operate with exactly one atom to measure the 
one-body heating. Second, we operate with ~ 3 trapped 
^^Rb atoms and use the result of the single-atom mea- 
surement to extract the two-body loss rate. To extract 
this rate we develop a Monte-Carlo simulation that we 
compare to the data. Our measurement indicates light- 
assisted collision rates that can reach remarkably large 
values (~ 10~^cm^.s~"'^), well above measured data found 
in the literature (by 1 — 2 orders of magnitude) for atoms 
held in magneto-optical traps (for a review, see [l|) or 
in larger dipole traps [3|, Ij] . Our maximal light-assisted 
collision rates are surprisingly close to the semi-classical 
Langevin limit. 



II. PRINCIPLE OF THE EXPERIMENT 

In our apparatus we operate an optical dipole trap at 
850 nm with micrometer size (waist 1 /im) [14| . which 
we load with cold atoms from a magneto-optical trap 
(MOT). Atoms enter the microscopic trap randomly, are 
trapped thanks to the cooling effect of the MOT beams, 
and are expelled from the trap due to one- or two-body 
processes. Depending on the local density of the MOT 
cloud around the dipole trap, we control the number of 
trapped atoms in steady state from one atom exactly 
(Nq = 1) to a few atoms on average ((iVo) ~ 3) [l7|. 

To study the light- assisted collisions, we switch off 
the MOT beams and then send the pulse of excitation 
light on the trapped atoms, initially prepared in the 
5S'i/27 F — 1 level. The excitation light consists of repum- 
ping light that transfers the atoms to the (55*1/2, F — 2) 
level, labeled S in Fig. [I][2]J, superimposed with light 




Figure 1: (a) Light-assisted atom loss by radiative escape. 
Atoms in the S state interact through the Van der Waals 
potential V{r) = —Ce/r^ while atoms in the S and P states 
interact through the dipole-dipole attractive potential V{r) = 
— Cs/r^. (b) Levels involved in the experiment. The excitation 
laser detuning S is measured with respect to the free space 
transition. 



nearly resonant with the {5Si/2,F — 2) ^ (5i^3/2,-F' ~ 
3) light-shifted transition, which excites the atoms into 
the (5P3/2, F' = 3) level, labeled P (see Fig. [lb). The ex- 
citation light consists of a pair of counter-propagating la- 
ser beams with orthogonal circular polarizations. For this 
experiment, we do not control the orientation of the ma- 
gnetic field, which has a magnitude smaller than 0.2 G. 

During the laser excitation, two atoms form a loosely 
bound pair with one atom in the S state and the other in 
the P state and interact through the long-range dipole- 
dipole attractive potential V{r) = —C^/r^ (here, r is the 
inter-atomic distance, C3 = 3?ir/4fc^, r/27r = 6 MHz is 
the linewidth of the P state, A = 27r/fc = 780 nm is the 
wavelength of the S ^ P transition), as represented in 
Fig.[T^. If the kinetic energy acquired by an atom pair be- 
fore it radiates back to the gound state exceeds the opti- 
cal dipole trap depth [/, it escapes the trap, thus leading 
to the loss of two atoms. This interaction-induced loss 
mechanism, known as radiative escape [15| [22] , co-exists 
with the standard one-body loss mechanism associated 
to the cycles of absorption and spontaneous emission of 
photons by individual atoms in the trap, which heat them 
out of the trap. Whether one or the other mechanism is 
dominant depends on the parameters of the experiment, 
namely the trap depth C/, the saturation s = ///sat (/ 
is the laser intensity and /gat = 1.6 mW/cm^) and the 
frequency detuning 5 of the excitation light with respect 
to the single atom transition in free space (see Fig. [TJd). 

For a given set of parameters we measure the number 
of trapped atoms that remain in the trap after the pulse 
of light has been sent. The number of atoms is measu- 
red by accumulating their fluorescence at 780 nm on an 
intensified CCD camera and com par ing it to the calibra- 
ted fluorescence of a single atom [l6|, |l7 1 . By varying the 
duration t of the pulse we obtain atom loss curves, from 
which we extract the one- and two-body loss rates as 
explained below. 



III. LIGHT-INDUCED LOSSES OF SINGLE 
ATOMS. 

To extract the one-body loss rate, we perform the loss 
experiment described above with one atom exactly in 
the trap (A^o — !)■ To do this, we adjust the loading 
rate of the trap to operate in the coUisional blockade 
regime [ij, [iSJ and trigger the loss experiment on the 
presence of a single atom in the trap. All parameters 
(trap depth, excitation light parameters) are otherwise 
unchanged with respect to the case {Nq) — 3 explored 
later in this paper. We obtain the survival probability 
of a single atom after the excitation process by repea- 
ting the experiment 200 times and measuring each time 
the presence or the absence of the atom in the trap after 
the experiment. Figure [5^1. shows examples of loss curves 
that illustrate the effect of the excitation light. For com- 
parison, the lifetime of the atom in the trap is 24 s in 
the absence of the excitation light and is limited by the 
residual background gas collisions. 

The effect of the excitation light is to heat the atom 
out of the trap as the duration of the excitation in- 
creases. This effect is quantitatively well explained by 
assuming that the temperature of the atom varies in time 
as T{t) = To + ai since the energy of the atom increases 
linearly with each absorption and spontaneous emission 
cycle. Here, Tq is the temperature at the beginning of the 
excitation pulse and a is the heating rate. Assuming a 
harmonic trap and a Boltzmann energy distribution |19| , 
the probability Pi{t) for a single atom to remain in the 
trap with depth t/ at a temperature T{t) is given by 



Piit) = l- 
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where r/(t) = U /k-BT{t) (fee is the Boltzmann constant). 
The temperature Tg being measured independently by a 
release-and- recapture technique [l^, we fit the data to 
the function Pi(t) with a being the only free parameter. 
The result obtained for a is shown in Fig. [5]d. As expec- 
ted, it reaches a maximum when the excitation light is 
nearly-resonant with the light-shifted S ^ P transition, 
i.e. 5 — U/h {h — 2-Kh is the Planck constant). 

The position and the shape of this resonance is confir- 
med by calculating the heating rate a = 2E^.R/k^, where 
R is the photon scattering rate and E^ is the recoil energy. 
We use a rate equation model to calculate the populations 
of the various Zeeman sub-levels from which we deduce 
the scattering rate R. We take into account the finite ini- 
tial temperature of the atom Tq ~ 200 /iK and its linear 
increase in time, which lead to random positions of the 
atom in the trap, and therefore to different light-shifts 
for the various Zeeman states. As seen in Fig. [SJ), this 
model reproduces the broadening of the resonance with 
respect to the natural line width r/27r = 6 MHz. In the 
following, we will use it to extrapolate the heating rate 
a to other values of Tq, as Tq varies when operating with 



o 



> 



3 
C/3 




"1 1 1 1 1 1 r 

10 20 30 40 50 60 70 
5/271 (MHz) 



Figure 2: (a) Survival probability Pi{t) of a single atom in 
the dipole trap, measured for various detunings of the ex- 
citation light 5/2n — {20; 40; 70} MHz (circles, squares and 
triangles, respectively). The trap depth is U = h x 36 MHz= 
fcB X 1.8 mK and the saturation parameter is s = 0.5. Solid 
lines show fits of the data to a heating model (see text). The 
initial temperature of the single atom is To — 200 ^K. Error 
bars are statistical, (b) Heating rate a (circles) deduced from 
the fits of the survival probability, versus 5/2-k (error bars are 
from the fit) ; the solid and dashed lines are solutions of the 
rate equation model described in the text, with no free para- 
meter, for To = 200 fiK and To = 300 ^K respectively. Data 
labeled 1, 2 and 3 are extracted from the loss curves shown 
in (a). 



IV. LIGHT-INDUCED LOSSES IN THE 
MULTI-ATOM CASE. 

A. Experimental observations 

We now turn to the case where a few atoms are loa- 
ded in the trap. Starting with (A^'o) — 3 atoms, we mea- 
sure the average number of atoms that remain in the 
trap after the excitation pulse has been sent [23| . Again 
the measurements are performed by averaging over se- 
veral hundreds of experiments. In the absence of exci- 
tation light the number of atoms remains constant on 
time scales large with respect to the pulse duration. In 
the presence of excitation light we observe losses that 
can be much faster than in the single atom case, depen- 
ding on the trap depth and the excitation light parame- 
ters. For example. Fig. [5^ compares loss curves taken for 
{Nq = l,ro = 200 ^K) and {{No) = 3.5, Tq = 300 ^K) 
initially, all other parameters being the same (s = 0.5, 
(5/27r = 20 MHz, U/h = 36 MHz). For this set of para- 
meters and (A'o) — 3.5, the number of atoms drops by a 
factor 2 in only 0.25 ms, at least one order of magnitude 
faster than in the single atom case. This rapid decrease is 
incompatible with the radiative heating rate a measured 
in the single atom case, taking into account the increase 
in To when we operate with a few atoms (see Sec. HI 
and Fig. ^jp). More generally, we observe this phenome- 
non for small values of the detuning, typically d ^ U/h 
(see Fig.IHJj). While a model involving only the radiative 
heating process does not reproduce the data, adding two- 
body losses to the model does, as shown in Fig. [3^. We 
thus attribute the observed excess losses in this regime 
to the leading two-body light-induced collisions, neglec- 
ting higher-body collisional processes. By contrast, for 
6 ^ U/h, loss curves overlap well in the few- and single- 
atom cases, indicating that the radiative heating process 
is identical in both cases and is the dominant loss me- 
chanism. In order to reveal the range of detunings where 
two-body losses dominate over one-body heating, we have 
represented in Fig. [SJd the inverse of the half-lifetime of 
the survival probability for both the single and the few- 
atom cases : for S/2Tr > 40 MHz « U/h, the two 
curves are nearly identical, indicating that the heating is 
dominant, while the two-body light- assisted losses domi- 
nate for smaller detunings. 



more than one atom initially. For instance, the dashed 
curve in Fig. [2}d shows how the resonance in a is shif- 
ted towards the low values of S when Tq increases from 
To = 200 fiK (as measured when Nq = 1) to Tq = 300 fiK 
(when (iVo) = 3). Qualitatively, this shift corresponds to 
the shift of the Boltzmann energy distribution towards 
the shallower parts of the trap. 



B. Model including radiative heating and two-body 

losses 



To extract the contribution of the two-body loss pro- 
cesses from our loss curves, we developed a Monte Carlo 
simulation to find the time-dependent number of atoms 
N{t). This approach is particularly appropriate in our si- 
tuation as the competing radiative heating process leads 
to a time-dependent one-body loss rate and because we 
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Figure 3: (a) Average number of atoms {N{t)) remaining in 
the dipole trap after excitation by a light pulse with dura- 
tion t. Squares : the initial number of atoms in the trap is 
(No) = 3.5 ; the solid line is a fit to a Monte-Carlo simula- 
tion (see text), with a heating rate a = 0.21 mK/ms and a 
two-body loss rate /3' = 1.2 (at.ms)"^. Circles : single-atom 
measurements rescaled to A'o = 3.5 ; the solid line is a fit to 
the radiative heating model. Error bars are statistical. Para- 
meters of the experiments : U/h = 36 MHz, S/2tt = 20 MHz, 
s — 0.5. (b) Inverse of the half-lifetime in the cases (A'^o) — 3.5 
(squares) and No = 1 (circles) versus the detuning of the ex- 
citation light. 



need to take into account the discreteness of the small 
atom number. We describe here the main lines of our 
simulation. 

At each time step dt we evaluate the infinitesimal pro- 
bability for a two-body loss event to have occurred within 
the N-atom ensemble between times t and t + dt, and we 
compare it to the probability of a one-body event to have 
occurred due to radiative heating during the same time 
interval dt. The first is denoted (ig'2-body and is related 
to the number of atom pairs at time t and the two-body 



loss constant (3' through [3| 

dq2.body^/3'N{t){N{t)-l)dt/2. 

The second is denoted d^i-body and is related to the num- 
ber of atoms N{t) and to the instantaneous one-body loss 
rate j{t) through 

rfgi-body = j{t)N(t)dt. 

Here, 7(f) = —Pi{t)/Pi{t), as obtained by a Taylor ex- 
pansion of Pi {t + dt) , where Pi (t) is given by Eq. [1] In 
practice, we calculate Pi{t) by using the heating rate 
a measured in the single- atom regime, corrected by the 
temperature Tq of the N-atom ensemble, which we mea- 
sure independently by a time-of-flight method. At each 
time step of the simulation, three channels are possible : 
(i) no loss occurs during dt : the probability associated to 
this channel is (1 — d9i-body)(l — rf92-body) ; (ii) a one-body 
loss takes place and the atom number decreases by one : 
the associated probability is d(7i-body(l — rf'?2-body) ; (iii) a 
two-body loss occurs and the atom number decreases by 
two : the associated probability is dg2-body(l — rfgi-body)- 
We pick up randomly one out of these three channels 
according to their associated probabilities, calculate the 
number of atoms at time t + dt, and then proceed to the 
next time step. By averaging over the initial atom num- 
ber distribution ^ 200 times we obtain a loss curve that 
simulates the actual measurements described above (see 
e.g. Fig. Oi). 

C. Light-assisted two-body loss rates 

The Monte Carlo simulation described in the previous 
section yields values of /3' that range from 0.02 (at.ms)"^ 
to 10(at.ms)~^, depending on the trapping and excita- 
tion parameters explored. In order to compare our results 
to theoretical models and to measurements reported el- 
sewhere in other trapping configurations, we calculate 
the normalized two-body loss rate f3 — (3' 2\/2V [l|, 
where V = C^^^T" )^^'^ is the volume occupied by the 
atoms assumed to be at thermal equilibrium at Tq. Here, 
Lj — (w^W||)^/'^ is the geometric average of the dipole 
trap oscillation frequencies uj± and wy, and m is the 
mass of an atom. For example, the data shown in Fig. [3^ 
are best fitted when /3' = 1.2 ± 0.5 (at.ms)"^. Using 
To = 300 /xK and the parameters of our set-up (wj^ = 
130 kHz and wy = 25 kHz), we obtain V = 0.7 finr^ and 
/3 = 2.4 ±1.1 X 10-9 cm3.s-i. 

We extracted in the same way the loss rate (3 for various 
values of the trap depth U and of the excitation detuning 
S and saturation s. Figure|4^ summarizes our results for a 
saturation parameter s = 0.5 when we scan the frequency 
of the excitation light across the trap depth. We observe 
a resonance in /? that is shifted to the red with respect to 
the frequency corresponding to the bottom of the trap, 
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Figure 4: Two-body loss rate versus (a) the light-shifted de- 
tuning S/2n — U/h, and (b) the saturation parameter s of the 
excitation light. Triangles, squares and lozenges correspond 
respectively to trap depths U/h — (18; 36; 54) MHz, initial 
peak atomic densities no = (2.7; 4.3; 5.1) x lO^'^ cm~'^ and ini- 
tial temperatures To = (200; 300; 400) /iK. The average initial 
atom number is (A^o) — 3 for all the data. Error bars are from 
the fits. The dashed lines are the Langevin limits Pl associa- 
ted to each set of parameters (f/, To) (see text). In (a), the 
saturation parameter is s = 0.5. In (b), the detuning of the 
excitation is S/2n — 10 MHz. The solid lines are fits of the 
data to a /3oo s/{l + s) model (with /3oo the only free parame- 
ter). The dotted lines indicate the asymptotical value /3cx) in 
each case. 



i.e. when the light-shifted detuning A = S — U/h ^ 0. 
On the blue-side of the resonance {5 > U/h), the two- 
body loss rate is suppressed, due to the excitation to a 
repulsive potential curve [ij, |20[ . 

We also observe that the peak value of /3 increases by 
more than an order of magnitude when the trap depth 
decreases only by a factor 3. Figure |4}d shows that /3 also 
increases as s/(l -|- s), in qualitative agreement with a 
simple model assuming a two-level system. For the largest 
saturation parameter investigated (s = 1.5), we find our 
largest value of the two-body rate constant f3 = 3.0 ± 



1.5 X IQ- 



D. Discussion 

Our measurements indicate that the light-assisted two- 
body loss rate can reach values remarkably larger than 
any reported measurements we could find, using three- 
dimensional excitation light and either ^^Rb or ^^Rb. For 
example, Kuppens et al. Q measure /3 ~ 10~^ cm'^.s"^ 
using a dipole trap with a waist of 26 /im and trap 
depths, temperatures and spatial densities comparable 
to ours. Kulatunga et al. [J] use a dipole trap with a 
waist of 5.6 /im size and measure two-body loss rates 
as large as /?' « 10^^ (at.s)^^. Estimating their volume 
at thermal equilibrium, we have found that this corres- 
ponds to a normalized loss rate (3 ^ 10~^^ cm'^.s"^. Only 



[a, |l8J have to assume large values of 
'^ to explain the loading of at most one 



Schlosser et al. 
13' - 1000 (at.s) 
atom in their sub- micrometer size dipole trap. Again esti- 
mating their one-atom thermal volume [SJ, ll9| , this yields 
/3~3x 10-9 cm3.s-i. 

The analysis presented in Sec. IV- C assumed the vo- 
lume V occupied by the atoms to be constant during 
the excitation pulse, and equal to the thermal volume for 
atoms at a constant temperature Tq. This assumption is 
actually not valid when the influence of the heating is lar- 
ger than or comparable to the two-body loss mechanism, 
i.e. when S <; U/h. However, neglecting the temperature 
increase during the light excitation leads to actually un- 
derestimate V and thus f3. For instance, we checked that 
for the highest values of /3 that we measured the tempe- 
rature increased by less than 15% during the excitation, 
leading to an underestimation of /3 by less than 25%, a 
difference within our error bars. 



E. Comparison to a semi-classical model 

Finally, we compare our largest measured light-assisted 
loss rate (i.e. ^oo = 3.1 ± 0.2 x 10"*^ cm^.s'^ in Fig.gjD) 
to the Langevin semi-classical limit f3i^ — ai^v 13| for the 
collision rate. Here, v = . l^hsTa jg ^j-^g average velocity 



of the atoms in the frame of the two-body center of mass. 
The Langevin cross-section ctl is obtained by summing 
the maximum cross-sections (2^ + l)47r?i /{mE) up to the 
maximum partial wave Zmax contributing to the collision, 
for a given collision energy E. In this approach, 
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We calculate Imax by imposing two conditions. First, as- 
suming the pair of atoms has been excited in the S + P 
potential (see Fig. [IJ, the kinetic energy E = ^ksTo of 
the two colliding atoms in the frame of their center of 



mass must be larger than the height of the centrifugal 
barrier to allow the collision to take place at short inter- 
atomic distance. This condition yields 
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as /max,i ^ 1 ('max,i — 50 typically). Second, the height 
of the centrifugal barrier in the S + S potential should be 
small enough to allow a pair of atoms with energy E to 
be excited at an interatomic distance shorter than their 
minimal approach distance. This yields 
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where the distance roxc actually depends on the light- 
shifted detuning A through r^xc — {—C^/hA)'^ . This se- 
cond condition yields 
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The maximal partial wave contributing to the cross- 
section is actually ^max = Min(?niax4,'max,2)- Figure H 
shows the Langevin limit set by these two conditions. 
For atoms at a temperature Tq — 200 ^K in a trap with 
U/h = 18 MHz and an excitation with 5/2tt = 10 MHz, 
/3l — 4.1 X 10^^ cm'^.s^^, close to our largest measured 
value /3oo = 3.1±0.2 x 10"*^ cm^.s"! (see Fig.Hb). Given 
the simplicity of the model, it is quite surprising that 
the light- assisted process studied here approaches this 
theoretical limit : for alkalies the two-body collision rate 
is predicted to be smaller than the Langevin limit by at 
least one order of magnitude [13 1. 



V. CONCLUSION 

In conclusion, using our ability to isolate one-body ra- 
diative heating from two-body losses, we have measured 
remarkably large two-body collision rates in a microme- 
ter size optical dipole trap in the presence of near re- 
sonant light. We have found that these large rates are 
close to the semi-classical Langevin limit. Given the com- 
plexity of the situation considered here, due to the near- 
resonant character of the light combined to the small size 
of the trapping potential that may affect the interaction 
between the atoms, it is quite remarkable that a simple 
semi-classical argument reproduces our largest measured 
value. It would be interesting to cross-check our findings 
using atoms in optical lattices, a situation where the sites 
also have a sub-micrometer size and where the number 
of atoms per site can be controlled precisely ll|, \l^ . 
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